AD/A-000  620 


INERTIAL  IMPACTION  EFFICIENCY  OF 
CYLINDRICAL  COLLECTORS  BY  DIGITAL 
TECHNIQUES  AND  EFFECTS  OF  PARTICLE 
SIZE  DISTRIBUTIONS 

Arthur  K.  Stu^mpfle,  et  a I 

Edge wood  Arsenal 

Aberdeen  Proving  Ground,  Maryland 


October  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


AD  A 0 00  6 


AD 


EDGEWOOD  ARSENAL  TECHNICAL  REPORT 
EC-TR-74050 

INERTIAL  IMPACTION  EFFICIENCY  OF  CYLINDRICAL  COLLECTORS 
BY  DIGITA  L TECHNIQUES  AND  EFFECTS 
OF  PARTICLE  SIZE  DISTRIBUTIONS 


Arthur  K.  Stuempfle 
Claude  R.  Pellegrino 


Chemical  Laboratory 


October  1 974 


I ;S  • A ft  Af  i 


u 


fc«0'^3uc*<l  hy 

fiATIOUAl  TECHNICAL 
(NFARVATIOM  ^FPVirF 

It  - O i* ft?  o*  Cr'-nm  rC< 


DEPARTMENT  OF  THE  ARMY 
Headquarters,  Edgewood  Arsenal 
Aberdeen  Proving  Ground,  Maryland  21010 


Approved  for  public  release;  distribution-unlimited. 


m 


- 6G0C20 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Pttm  Entered)  , 

REPORT  DOCUMENTATION  PAGE  LfoI^cohIl^gform 

r REPORT  NUMBER  |2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

F.C-TR-74050 


4.  TITLE  (end  Subtitle)  5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

INERTIAL  IMPACTION  EFFICIENCY  OF  CYLINDRICAL  Technical  Report 
, COLLECTORS  BY  DIGITAL  TECHNIQUES  AND  EFFECTS  February-May  1974 
I OF  PARTICLE  SIZE  DISTRIBUTIONS  6.  performing  org.  report  number 


7.  AUTHORf*; 

Arthur  K.  Stuempfle 
Claude  R.  Pellegrino 

9.  performing  organization  name  and  address 
Commander,  Edgewood  Arsenal 
Attn:  SAREA-CL-POA 
Aberdeen  Proving  Ground,  Maryland  21010 

1 1.  CONTROLLING  OFFICE  name  ano  address 

Commander,  Edgewood  Arsenal 
Attn:  SAREA-TS-R 

Aberdeen  Proving  Ground.  Maryland  21010 

14.  MONITORING  AGENCY  NAME  a AODRESSf//  different  trom  Controlling  Office) 


| B.  CONTRACT  OR  GRANT  NUMBERfaJ 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS 


1W1621 16A08402 
__________  — - — - 

October  1974 

13.  NUMBER  OF  PAGES 

4L 

IS.  SECURITY  CLASS,  (ol  thle  report) 

UNCLASSIFIED 

ISa.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 

NA 


16.  DISTRIBUTION  STATEMENT  ol  this  Report < 


Approved  for  public  release:  distribution  unlimited. 


I 17.  DISTRIBUTION  STATEMENT  *1  the  abstract  entered  In  Mock  20.  li  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


Chemical  test  and  assessment  technology 


1 19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


Droplet 

Aerosol 

Sampling 

Particles 

Impaction 


Cylinders 

Collectors 

Deposition 

Efficiency 

Impingement 


Interception 
Particle  impaction 
Inertial  impaction 
Inertial  parameter 


Sampling  efficiency 
Impaction  efficiency 
Collection  efficiency 
Deposition  efficiency 


Isokinetic  sampling  Weibull  distribution 


20.  ABSTRACT  (Continue  on  reveree  title  II  necessity  end  identity  by  block  number) 

The  theory  of  inertial  impaction  ol  particles  on  cylinders  has  been  analyzed  and  used  to  develop 
standard  inertial  impaction  efficiency  curves.  UiiKpic  exponential  functions  have  been  generated  by  a 
digital  computer  that  accurately  fit  the  inertial  impaction  theory  with  a maximum  relative  error  of 
less  than  1%.  Interpolate  routines  have  been  adapted  for  the  computer  program  to  obtain  inertial 
impaction  efficiency  predictions  lor  all  incilt.il  parameter  and  velocity  field  scaling  parameter  values 
in  the  range  of  0.1 3 < K < 300  and  0 < <t>  < 10.000.  The  complete  computer  program  with  examples 
and  solutions  of  test  cases  are  presented  in  the  appendix. 

(Continued) 


DD  ,^nrm73  1473 


EDITION  OF  I NOV  65  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wien  Dora  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  FAOEOnxn  Dim  Kntmnd) 


19.  KEYWORDS 

Drop  size  distribution 
Digital  computet  program 

20.  ABSTRACT 

The  Weibul!  distribution  function  has  been  adapted  to  analyze  the  effects  of  particle  size 
distributions  on  the  impaction  efficiency  of  cylinders.  Use  of  the  mass  median  diameter  to 
characterize  a particle  size  distribution  is  unsatisfactory  for  predicting  the  impaction  efficiency 
from  heterogeneous  aerosols.  Use  of  the  particle  size  distribution  to  predict  efficiency  of 
impaction  yields  good  agreement  between  theory  and  experiment. 


■ UNCLASSIFIED 

SECURITY  CLASSIFICATION  >F  YHIS  PAGEfWitfrt  i>*»*  hntetod) 


tail 


0 > ' 


PREFACE 


The  work  described  in  this  report  was  authorized  under  Tasl  1W1621 16A08402,  Chemical  Test  and 
Assessment  Technology.  This  work  was  started  in  February  and  completed  in  May  1974. 

Reproduction  of  this  document  in  whole  or  in  part  is  prohibited  except  with  permission  of  the 
Commander,  Edgewood  Arsenal,  Attn.  SAREA-TS-R,  Aberdeen  Proving  Ground.  Maryland  21010.  however,  DDC 
and  the  National  Technical  Information  Service  are  authorized  to  reproduce  this  document  for  United  States 
Government  purposes. 


r 


i 


CONTENTS 


I.  INTRODUCTION 

II.  BACKGROUND 

III.  METHODS  AND  RESULTS 

IV.  DISCUSSION 

V.  CONCLUSIONS 

LITERATURE  CITED 

APPENDIX,  Digital  Computer  Program  for  Inertial  Impaction  Efficiency  Computations  . . 
DISTRIBUTION  LIST 


Page 

5 

5 

12 

17 

19 

21 

22 

39 


INERTIAL  IMPACTION  EFFICIENCY  OF  CYLINDRICAL  COLLECTORS  BY  DIGITAL  TECHNIQUES 
AND  EFFECTS  OF  PARTICLE  SIZE  DISTRIBUTIONS 


I.  INTRODUCTION. 


The  deposition  of  particles  of  intermediate  size  ( 10  to  200  pm  diameter)  on  collectors  in  a flow  stream  is 
principally  based  on  the  inertia]  impaction  mechanism.  Incttial  impaction  of  panicles  on  collectors  is  of  interest  in 
diversified  areas  such  as  crop  dusting,  mosquito  spraying,  an  pollution  contiol,  aiiciafi  icing,  oaiticuiate  buildup  on 
heat  exchangers,  and  in  practically  any  circumstance  where  matter  in  particulate  foim  is  removed  from  a tiansport 
fluid.  The  efficiency  with  which  the  particles  are  removed  from  the  flow  stream  is  a function  of  the  particle  si/.c, 
collector  size,  and  flow  field  conditions.  For  a stationary  cylinder  in  a moving  airstrcam,  the  impaction  efficiency  is 
defined  as  the  ratio  of  the  cross-sectional  area  of  the  upstream  envelope  containing  the  trajectories  of  the  particles 
which  intersect  the  collector  surface  to  the  cross-sectional  area  of  the  cylindrical  collector  norma!  to  the  direction  of 
flow.  Estimates  of  the  impaction  efficiency  for  a given  set  of  conditions  can  be  found  by  computing  the  trajectories 
of  the  particles  that  challenge  the  collector  or  by  graphically  determining  the  efficiency  from  curves  constructed 
from  the  pomt-by -point  trajectory  calculations.  These  methods  »ie  adequate  when  only  a few  data  points  arc  of 
interest  but  the  technique  becomes  tedious  when  estimating  efficiencies  for  a variety  of  possible  impaction 
conditions.  In  addition,  the  aerosol  that  challenges  the  collector  generally  consists  of  a wide  range  of  sizes  depending 
on  the  method  used  for  particle  generation.  It  is  recognized  that,  if  a single  parameter,  such  as  the  mass  median 
diameter  (MMD),  is  used  to  characterize  the  particle  size  distribution  for  computational  purposes,  the  expected 
theoretical  impaction  efficiency  can  be  grossly  different  from  the  experimental  efficiency.  These  differences  can  be 
attributed  in  part  to  the  turbulent  nature  of  the  flow  field  and  to  the  uonliueai  character  of  the  inertial  impaction 
efficiency  curves.  The  objectives  of  this  study  have  been  to  develop  a computer  program  that  rapidly  calculates  the 
inertial  impaction  efficiency  of  cylindrical  collectors  for  any  given  set  of  laminar  flow  conditions  and  to  devise  a 
technique  to  assess  the  effects  of  a particular  particle  size  distribution  on  impaction  efficiency . 

II.  BACKGROUND. 


The  theory  of  inertial  impaction  of  spherical  particles  on  cylindrical  collectors  in  an  ideal  flow  field  is 
based  on  a numerical  solution  of  the  equations  of  motion  of  the  particles  undergoing  tiansport  uiound  the  bluff 
body.  Development  of  the  theory  and  associated  digital  computer  techniques  lias  been  described  in  a previous  study  ' 
and  provides  the  basis  for  this  report.  Langmuir  and  Blodgett2  proposed  the  following  dimensionless  form  of  the 
equations  of  motion  for  a particle  in  Cartesian  coordinates. 
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vx.  v = particle  velocity  components  normalized  by  the  tree-stream  velocity 
Ujj,  u = airstrcam  velocity  components  normalized  by  the  Iree-stieam  velocity 

U = free-stream  velocity  at  an  infinite  distance  from  the  cylmdei  surface 


( c 


D 


Re  = 


K = 


P 

dP_ 

pd2u 

18pR 


P 


drag  coefficient  for  spherical  particles  in  fluid 

particle  Reynolds  number  with  respect  to  local  relative  velocity 
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R = cylinder  radius 

r = — t = time  scale 
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The  inertial  parameter,  K,  is  a measure  of  the  inertia  of  the  particle  and  relates  to  the  magnitude  of  the  external 
force  required  to  cause  a change  in  its  direction  of  motion.  If  Stokes’  law  of  resistance  is  assumed,  the  inertial 
parameter  represents  the  ratio  of  the  “stopping"  distance  of  a particle  projected  with  velocity  U into  still  air  to 
the  radius  of  the  cylinder.  A second  dimensionless  parameter  proposed  by  Langmuir  and  Blodgett  is  the  velocity 
field  scaling  parameter,  0,  originally  defined  ir.  terms  of  the  free-stream  particle  Reynolds  number.  The  <•> 
parameter  is  used  to  calculate  the  magnitude  of  the  instantaneous  Reynolds  number  of  the  particle  at  any  point 
in  the  flow  field  and  can  be  related  to  the  Reynolds  number  of  the  cylindrical  collector  as  follows: 
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where 

Va° 

ReQ  = — ^ — = free-stream  Reynolds  number  of  the  particle 

Rec  = free-stream  Reynolds  number  of  the  cylindrical  collector 

The  magnitude  of  the  <j>  parameter  at  any  point  in  the  particle  path  is  a measure  of  the  deviation  from  Stokes’ 
law  due  to  the  forces  acting  on  the  particle. 


The  starting  point  for  a mathematical  description  of  the  fluid  flow  field  around  the  collector  is  the 
Navier-Stokes  equation.  A solution  to  the  general  Navier-Stokes  equation  is  possible  after  a number  of  limiting 
assumptions  are  made  with  respect  to  the  fluid  and  flow  conditions.  One  form  of  solution  can  be  found  if  the 
fluid  is  assumed  to  be  ideal;  i.e.,  constant  density,  irrotational,  and  without  viscosity.  The  flow  pattern  around 
the  cylinder  under  steady-state  conditions  is  dependent  on  the  Reynolds  number  of  the  collector.  Where  the 
Reynolds  number  of  the  cylinder  is  one  thousand  or  greater  (0  > 10  in  air),  the  flow  can  be  considered  ideal  in 
the  absence  of  turbulence  and  the  flow  field  is  adequately  described  by  potential  theory  for  an  incompressible 
fluid.  The  streamlines  for  flow  outside  the  cylinder  can  be  obtained  and  the  airstream  velocity  components  can 
be  written  in  dimensionless  terms  as  simple  functions  of  the  reduced  position  coordinates’;  namely. 
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Figure  1.  Coordinate  System  for  Cylinder  in  Potential  Flow  Field 


where  the  Cartesian  coordinates  have  been  normalized  with  respect  to  the  cylinder  radius. 


A numerical  solution  of  the  differential  equations  of  motion  given  above  requires  the  instantaneous  drag 
force  on  the  particle  to  be  estimated  as  the  trajectory  is  developed.  The  steady-state  drag  coefficient  of  the  sphere  is 
a function  of  the  Reynolds  number,  and  the  reliable  experimental  da.a  as  tabulated  by  Fuchs3  have  been 
interpolated  to  obtain  these  estimations.  The  probable  error  in  the  drag  coefficient  for  solid  spheres  is  claimed  by 
Fuchs  to  be  less  than  4%  over  the  Reynolds  number  range  of  0.01  to  500.  These  data  are  nearly  identical  to  the 
experimental  data  cited  by  Schlichting4  as  adopted  and  reported  by  Hussein  and  Tabakoff.5  Over  the  Reynolds 
number  range  of  0.10  to  7.0,  the  earlier  approximations  used  by  Langmuir  and  Blodgett  differ  from  the  current  drag 
coefficient  data. 

The  Cartesian  coordinate  system  used  in  calculating  the  trajectory  of  particles  around  the  cylindrical 
collector  is  identical  to  that  of  Brun,  Lewis,  Perkins,  and  Serafini®  and  is  shown  in  figure  1.  The  motion  of  the 
particles  is  in  a plane  perpendicular  to  the  cylinder  axis  which  is  the  origin  of  the  coordinate  system.  Theoretically, 
the  maximum  angle  of  impinginent,  is  the  angle  beyond  which  no  deposition  occurs  by  the  inertial  mechanism. 

Several  simplifying  assumptions  have  been  made  in  the  derivations  and  in  computing  the  trajectory  of  a 
particle.  These  assumptions  include  the  following: 

1.  At  an  infinite  distance  upstream  of  the  cylinder,  the  particles  have  horizontal  and  vertical  velocity 
components  equal  to  the  free-stream  air. 

f 2 The  particles  are  spherical,  monometric  (uniform  size),  and  monodisperse  (single  particles),  and  they 

do  not  evaporate  or  deform. 

3.  Gravitational,  electrostatic,  and  any  other  external  forces  are  negligible. 

4.  The  particle  radius  is  negligible  with  respect  to  the  cylinder  radius.  (For  interception  effect 
considerations,  see  previous  publication.7) 

\ 5.  The  boundary  layer  about  the  cylinder  surface  does  not  affect  the  particle  trajectory. 

\ 6.  The  airflow  around  the  cylinder  is  described  as  ideal  and  without  circulation  and  is  unaffected  by  the 

presence  of  the  particles. 

! 

; 7 The  instantaneous  drag  force  coefficient  for  the  particle  is  given  by  the  steady-state  data  and  is  not 

subjected  to  acceleration  effects. 

8.  All  particles  that  strike  the  collector  adhere  to  it. 

\ 

Results  of  the  digital  computer  trajectory  calculations  to  determine  the  efficiency  of  impaction  of 

‘ spherical  particles  on  cylindrical  collectors  in  an  ideal  flow  field  as  a function  of  the  K and  <p  scaling  parameters  are 

shown  in  figure  2.  The  impaction  data  for  inertial  parameter  values  less  than  one  (K  < 1)  arc  displayed  in  figure  3. 

i In  addition  to  the  impaction  efficiency,  the  maximum-angle-of-impingemcnt  data  are  plotted  in  figure  4.  The  theory 

has  been  experimentally  verified  under  laminar  flow  conditions  for  a wide  range  of  inertiai  parameter  and  velocity 

. field  scaling  parameter  values.8 * *'1 1 In  general,  however,  the  theory  is  inappropriate  for  accurately  predicting  the 

impaction  efficiency  of  particles  possessing  small  inertial  parameter  values  on  collectors  in  fluid  flows  that  exhibit 
turbulence  intensity  levels  above  7.5 % or  when  the  Reynolds  number  of  the  collector  is  less  than  approximately  one 
thousand  ( <p  ^ 10  in  air).11  Further,  the  collection  efficiency  of  a cylinder  that  is  challenged  by  a distribution  of 
particle  sizes  is  not  accurately  predicted  by  a single  particle  parameter  such  as  the  mass  median  diameter.9  13  For 
many  practical  circumstances,  though,  the  theory  is  quite  sufficient,  and  this  study  has  been  undertaken  so  that 
impaction  efficiency  predictions  can  be  raplJly  made  by  computer  with  minimum  input  data  requirements  or 
computer  programming  background. 
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III.  METHODS  AND  RESULTS. 

Results  of  the  point-by-point  trajectory  calculations  define  the  inertial  impaction  theory  as  shown  in 
figures  2 and  3 and  provide  the  impaction  efficiency  as  a function  of  the  K and  $ scaling  parameters.  These  data  serve 
as  the  standards  for  which  unique  exponential  functions  have  been  obtained  to  fit  the  curves.  The  analytical 
expressions  for  the  standard  $ curves  can  then  be  interpolated  to  determine  the  impaction  efficiency  for  any 
nonstandard  K and  <j>  value.  A digital  computer  program  has  been  developed  to  perform  these  calculations  and  only  a 
minimum  of  input  data  is  necessary  for  operation. 

The  standard  velocity  field  scaling  parameter  curves  selected  for  fit  include  $ =0;  1;  2.2;  5;  10;  22;  50; 
100;  200;  500;  1,000;  2,210;  5,000;  and  10,000.  Each  curve  has  been  subdivided  into  four  inertial  parameter  ranges 
and  each  portion,  respectively,  has  been  fitted  by  the  same  form  of  the  unique  exponential  functions.  Any 
differences  that  occur  between  the  computed  efficient  along  a standard  £ curve  and  the  point-by-point  impaction 
efficiency  data  are  in  the  third  or  higher  significant  digit  position  and  the  maximum  relative  error  observed  has  been 
less  than  1%. 


Lagrange’s  interpolation  formula  has  been  used  to  compute  impaction  efficiency  values  for  velocity  field 
scaling  parameters  enclosed  within  the  aforementioned  standard  tj>  values  that  are  equal  to  or  greater  than  one.  The 
formulation  used  in  this  study  has  been  defined  as 


where 
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and 

0.13  < K < 64 


0j  = 1 ; 2.2;  5;  10; . . . 5,000;  10,000  - 

i,  j = 1,  2,  ...,  13  . i 

M = 1,  2 13 

• « 
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■| 


If 64  < K < 300, 

= 1;5;22;  100;  500;  2,210;  10,000 

i,  j - 1,2, ...,  7 
M = 7 


where  K = inertial  parameter  value  of  interest;  0.13  < K < 300. 

One  advantage  of  Lagrange’s  interpolation  formula  over  other  interpolative  routines  is  that  equidistance  values  of  the 
independent  variable  are  not  necessary  to  obtain  a reliable  value  of  the  dependent  variable.  For  4>  values  between  0 
to  1 , where  the  potential  flow  assumption  may  not  be  valid,  a b'near  interpolation  has  been  employed. 

Comparisons  have  been  made  between  impaction  efficiency  values  obtained  by  the  interpolative  method, 
and  the  point-by-point  trajectory  technique  and  relative  errors  are  less  than  1%  when  the  absolute  impaction 
efficiency  is  greater  than  5%.  For  impaction  efficiencies  less  than  5%,  the  largest  relative  error  observed  has  been 
within  10%  of  the  point-by-point  trajectory  calculated  efficiency  value  and  is  due  principally  to  round-off 
(four-place  accuracy).  The  trajectory  calculations  had  been  terminated  at  the  inertial  parameter  value  of  64. 
However,  the  interpolative  routine  extends  this  range  to  K values  of  300  where  the  impaction  efficiency  curves 
asymptotically  approach  the  theoretical  limiting  efficiency  of  100%  for  point-mass  particles.  Consequently,  for  all 
practical  purposes,  the  interpolative  routine  provides  an  accurate  prediction  of  the  impaction  efficiency  of  spherical 
particles  on  cylindrical  collectors  over  an  extensive  range  of  inertial  parameter  values.  A complete  and  detailed 
description  of  the  digital  computer  program  developed  in  this  study  is  contained  in  the  appendix. 

Experimental  verification  of  the  inertial  impaction  theory  has  been  accomplished  over  a wide  range  of 
parameter  values.  The  theory  turns  out  to  be  remarkably  accurate  for  the  monometric  (single-size)  particles  used  in 
the  experimental  efforts  especially  in  view  of  the  assumptions  that  have  been  made  regarding  fluid  properties  and 
flow  conditions. 

In  practical  applications,  however,  the  aerosol  that  impacts  on  a collector  will  seldom  consist  of  particles 
of  uniform  size.  Tire  size  distribution  of  the  aerosol  is  a function  of  the  methods  and  techniques  used  in  generation 
of  the  particles.  A single  empirical  or  theoretical  equation  does  not  exist  that  can  universally  predict  the  particle  size 
distribution  resulting  from  the  dispersion  of  liquids  by  spray  nozzles,  hot  and  cold  gas  atomizers,  and  explosive 
disseminators  or  the  dispersion  of  solids  from  mechanical  dispersers. 

A vast  literature  has  been  developed  by  various  authors  in  their  attempts  to  obtain  mathematical 
distribution  functions  that  characterize  experimentally  generated  particle  size  distributions.  Familiar  examples 
include  the  Nukiyama-Tanasawa  equation  for  drop  sizes  in  sprays  generated  by  air  atomization,  the  Rosin-Rammler 
equation  for  pulverized  coal,  Roller’s  formula  for  powder  materials,  the  normal  distribution  function  for 
symmetrically  distributed  particles  of  narrow  size  range  such  as  plant  spores,  and  the  log-normal  distribution  for  a 
large  number  of  condensation,  natural,  and  mechanically  generated  aerosols.  However,  the  intent  of  this  study  is  not 
to  generalize  on  particle  size  distributions  or  delve  into  their  specific  merits  and  methods  but  rather  to  examine  the 
necessity  of  accounting  for  the  particle  size  distribution  when  predicting  the  impaction  efficiency  of  a particular 
cylindrical  collector. 
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A simple  technique  used  to  obtain  a size  distribution  has  been  to  microscopically  observe  and  measure  a 
representative  sample  of  an  aerosol  population.  The  number  of  particles  that  lie  between  radius  r and  r + dr  can  be 
found  as  a function  of  the  radius.  The  fraction  of  the  total  number  of  particles  that  lie  in  an  interval  is  given  as 
an  - f£«)dr 


where 


CO 


i(r)dr  = 1. 


The  relative  frequency  distribution  curve  is  found  by  plotting  the  resulting  data  representing  dn,  and  the 
various  distribution  functions  mentioned  above  can  provide  a mathematical  expression  of  the  function,  f(r). 

Distribution  functions  can  be  mathematically  expressed  and  used  in  a number  of  ways  depending  on  the 
problem  requirements.  In  some  instances,  the  cumulative  fraction  of  particles  having  radii  greater  or  less  than  some 
radius  is  convenient.  Thus,  respectively. 


F(r)  = 


f(r)dr  or 


r 

f(r)dr 


In  addition  to  the  number  distribution,  the  mass  distribution  function  can  be  determined  and  is  useful  in  many 
practical  applications.  The  mass  fraction  of  particles  having  radii  between  x and  x + dx  is  written  as  df  = f(x)dx 


and  I 
Jo 


f(x)dx  = 1.  The  mass  distribution  function  is  related  to  the  number  distribution  function  as 

f(x)  = to  m,ffr) 


where 


mf  = mass  of  particle  with  radius  r 
co  = proportionality  constant 


With  respect  to  impaction  of  particles  on  collectors,  one  is  usually  interested  in  the  mass  distribution 
function  because  the  dose  acquired  by  the  collector  is  dependent  on  the  aerosol  mass  that  is  deposited  over  the 
sampling  period.  In  order  to  estimate  the  mass  collected,  an  average  particle  size  is  generally  computed  and  applied 
to  the  inertial  impaction  theory.  The  objectives  of  using  an  average  diameter  or  other  measure  of  central  tendency 
are  to  provide  a single  number  that  will  simply  desvrite  the  behavior  of  the  entire  aerosol  population  and  to 
eliminate  an  extraordinary  number  of  tedious  calculations.  The  problem  is  that  the  inertial  parameter  value 
computation  makes  use  of  the  square  of  the  particle  diameter  and,  further,  the  impaction  curves  are  a nonlinear 
function  of  the  inertial  parameter.  Consequently,  an  accurate  estimation  of  impaction  efficiency  or  mass  deposit  can 
only  be  expected  when  a single  particle  size  estimator  is  used  if  the  particle  size  distribution  is  over  a very  narrow 
range. 


In  order  to  demonstrate  the  effect  that  a particle  size  distribution  has  on  the  impaction  efficiency, 
experimental  data  are  required  for  comparison  purposes.  Most  of  the  experimental  work  that  has  been  performed  in 
the  past  to  verify  the  inertial  impaction  theory  has  been  conducted  with  monometric  (single-size)  aerosols.  One 
notable  exception  is  that  of  Landahl  and  Herrmann.12  Several  authors  refer  to  these  data  but  they  consider  the 
results  to  be  of  limited  value  due  to  the  use  of  heterogeneous  aerosols,  especially  at  small  values  of  the  inertial 
parameter. 
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Landahl  and  Hermann  studied  the  deposition  of  aerosols  on  wires  and  cylinders  with  particles  generated 
by  an  impinger-type  atomizer  and  by  a “modified  Brink’s  nozzle.”  The  size  distributions  v/ere  obtained  by  cascade 
impactorj,  and  microscopic  observations  and  representative  distributions  of  the  clouds  produced  are  illustrated  in 
the  following  table: 


Diameter 

Percent* 

Small  particle 

Medium  particle 

Large  particle 

cloud 

cloud 

cloud 

pm 

pm 

pm 

5 

(0.8) 

1.8 

5 

10 

1.0 

2.9 

8 

50 

4.0 

13 

28 

90 

12 

40 

58 

98 

(20) 

55 

75 

♦Percent  of  mass  less  than  stated  diameter. 


These  data  arc  not  easily  fitted  or  well  represented  by  the  aforementioned  distribution  functions. 
However,  a general  statistical  distribution  function  that  has  found  wide  applicability  recently  is  the  Weibull 
distribution  function.13"15  This  distribution  function  has  been  successfully  used  in  processes  involving  limits 
and  maxima/ininima  problems  that  include  lifetime  and  failure  rate  distributions  of  electrical  systems.  It  has  also 
been  used  for  describing  bounded  particle  size  distributions. 

The  cumulative  distribution  function  of  the  Weibull  distribution  has  been  applied  to  the  Landahl 
and  Herrmann  data  and  has  been  expressed  in  the  following  form: 


where 


F(x)  = 1 * exp 


(<>) 


F(x)  = cumulative  mass  fraction 

x = particle  diameter  (pm) 

7 = minimum  size;  location  parameter  (pm) 

j?  = characteristic  size;  scale  parameter  (pm) 

0 = shape  parameter 


The  triai-and-error  method  of  estimating  the  Weibull  parameters  has  been  described  by  Nelson15  and 
a digital  computer  program  has  been  developed16  to  estimate  the  parameters  of  the  distribution  function  for  the 


15 


\ 

l 

I 


small,  medium,  and  large  particle  clouds.  The  parameter  approximations  obtained  by  the  principle  of  maximum 
likelihood  method  turn  out  to  be  as  follows: 


Parameter 

Small  cloud 

Medium  cloud 

Large  cloud 

V 

5.3053 

17.9678 

34.4356 

P 

0.9532 

1.1183 

1.7461 

7 

0.5919 

0.5269 

-1.7375 

Evaluating  the  cumulative  mass  function  for  the  experimental  particle  diameters  to  assess  the  Weibull  distribution 
fit  yields  the  following  comparisons: 


Experimental 

F(x) 

Small  cloud 

Medium  cloud 

Large  cloud 

X 

F(x)  calculated 

X 

F(x)  calculated 

X 

F(x)  calculated 

% 

/am 

% 

/am 

% 

/am 

% 

5 

0.8 

5 

1.8 

5 

5 

5 

10 

1.0 

9 

2.9 

10 

8 

10 

50 

4.0 

52 

13 

50 

28 

51 

90 

12 

90 

40 

92 

58 

91 

98 

20 

98 

55 

97 

75 

98 

As  observed,  an  adequate  fit  is  obtained  for  ali  clouds  considering  the  limited  available  data  and  the 
wide  size  distributions  that  are  involved.  It  must  be  emphasized  again  that  the  objective  is  not  to  find  the 
“correct”  distribution  function  to  describe  a particular  particle  size  distribution  but  rather  to  approximate  the 
function  as  simply  as  possible  and  to  test  whether  the  resultant  distribution  significantly  affects  the  impaction 
efficiency  predictions. 

Given  the  above  distribution  functions  and  approximate  values  for  the  Weibull  parameters,  one  can 
solve  equation  9 for  the  particle  size  as  a function  of  the  cumulative  mass  fraction.  F(x);  namely, 

X ■ r + On-)  [m^  )]'"  (10) 


Subsequently,  by  selecting  appropriate  values  for  F(x),  the  average  particle  size  associated  with  a 
given  particle  mass  interval  can  be  determined  by  difference.  For  example,  if  30%  of  the  cloud  mass  has  particle 
sizes  below  7.4  /am  and  31%  of  the  cloud  mass  has  particles  less  than  7.6  /im  in  diameter,  then  1%  of  the  cloud 
mass  is  represented  by  particles  with  an  average  diameter  of  approximately  7.5  /am.  The  impaction  efficiency  for 
each  size  interval  can  then  be  computed  and  an  average  efficiency  can  be  obtained  for  the  entire  distribution. 
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A computer  program  has  been  written  to  perform  the  tedious  computations  for  each  1%  mass 
interval,  and  the  average  distribution  impaction  efficiency  is  found  based  on  96"  of  the  cloud  mass.  In  addition, 
the  program  provides  the  impaction  efficiency  based  stnctly  on  the  mass  median  diameter  of  the  cloud.  The  only 
supplementary  imput  requirements  necessary  for  use  of  the  distribution  function  option  to  the  computer 
program  are  the  Weibull  parameters  jj,  0,  7. 

The  sample  distribution  functions  have  been  used  to  compute  average  impaction  efficiencies  for 
comparison  with  the  experimentally  obtained  efficiencies  of  Landahl  and  Herrmann.  The  resuits  of  these 
computations  are  included  in  table  1.  The  velocity  field  scaling  parameter  values  for  the  Landahl  and  Herrmann 
tests  are  very  low  and  theoretically  the  potential  flow  assumptions  regarding  the  fluid  flow  field  should  not 
apply.  Further,  Landahl  and  Herrmann  state  that  the  airflow  was  turbulent,  but  the  degree  ot  turbulence  is  not 
indicated.  In  addition,  the  illustrative  particle  size  distributions  provided  by  Landahl  and  Herrmann  do  not 
correspond  with  their  reported  efficiency  data  except  for  the  13-pm-MMD  case.  .Nevertheless,  examination  of 
table  1 clearly  shows  that  the  theoretical  impaction  efficiencies  for  the  distribute'!:  functions  more  closely 
predict  the  observed  impaction  effii.ieni.ics  than  the  mass  median  diametci  computed  efficiencies  for  almost  al! 
circumstances.  This  is  especially  tru.  for  the  13-pm-MMD  case  where  the  assumed  particle  size  distribution 
apparently  matches  the  experimental  mass  median  diameter  and  the  theoretical  efficiencies  turn  out  to  be 
considerably  more  accurate  than  the  mass  median  diameter  prediction.  Note  that  the  small  and  large  <loud 
particle  size  distributions  have  mass  median  diameters  that  exceed  the  experimental  mass  median  diameters  and. 
therefore,  tend  to  overestimate  the  theoretical  tffi.iencies  from  the  distribution  function.  Thus,  to  summarize  the 
results,  when  the  particle  size  distribution  is  used  to  compute  the  average  impaction  efficiency  for  the  cloud,  the 
prediction  corresponds  much  more  closely  to  the  observed  efficiency  than  that  of  a simple  mass  median  diameter 
prediction.  The  predictions  arc  relatively  accurate  when  the  particle  size  distribution  is  known  as  evidenced  by 
the  13-pm-MMD  test  case. 

IV.  DISCUSSION. 

The  important  conclusion  that  can  be  drawn  from  the  theoretical  and  experimental  data  presented  in 
table  1 is  that  use  of  the  mass  median  diameter  to  describe  a particle  size  distribution  will  seldom  result  in  an 
accurate  prediction  of  the  impaction  efficiency  for  the  cloud.  In  general,  at  large  values  of  the  inertial  parameter. 
K,  the  mass  median  diameter  overestimates  the  impaction  efficiency  of  the  distribution,  whereas,  for  small  values 
of  the  parameter  where  the  theoretical  efficiency  for  the  mass  median  diameter  is  less  han  1 7i,  the  efficiency 
can  be  underestimated  by  orders  of  magnitude.  These  differences  arc  due  to  the  nonlinear  character  of  the 
inertial  impaction. theory.  However,  when  tire  particle  size  distribution  is  considered  in  the  predictions,  the 
agreement  between  theory  and  experiment  is  much  better. 

Calculations  were  performed  using  the  typical  particle  size  distributions  previously  given  to 
determine  the  range  of  inertial  parameter  values  over  which  the  mass  median  diameter  predictions  were  within 
2:10%  of  the  average  impaction  efficiencies  for  the  distributions.  In  general,  if  the  mass  median  diameter  yields  a 
K value  between  1 and  2,  for  all  $ values,  the  effect  of  the  distribution  on  efficiency  need  not  be  considered. 
When  the  mertial  parameter  value  based  on  the  mass  median  diameter  lies  outside  this  range,  the  average 
distribution  efficiency  will  differ  from  the  mass  median  diameter  predicted  efficiency  by  more  than  10".  To 
accurately  predict  the  impaction  efficiency  for  heterogeneous  aerosols  when  the  K value  for  the  mass  median 
diameter  is  less  than  one,  the  particle  size  distribution  must  be  taken  into  account.  Unfortunately,  data  arc  not 
available  to  determine  the  degree  of  monometncity  of  the  distribution  that  is  required  in  order  that  the  mass 
median  diameter  or  other  average  particle  diameter  will  adequately  represent  the  population  over  the  entire 
inertial  impaction  region  of  interest. 
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It  is  also  interesting  that  the  predicted  efficiencies  arc  close  to  die  experimentally  measured 
efficiencies  even  though  the  velocity  scaling  parameter  values  arc  much  less  than  10  in  some  cases.  In  general,  the 
lower  <p  values  would  lead  to  reduced  impaction  efficiencies  because  the  flow  field  would  tend  to  exhibit  viscous 
flow  properties.  On  the  other  hand,  no  account  of  interception  effects  has  been  made  for  the  interpolation 
routine,  and  interception  could  occur  since  the  particle  to  collector  diameter  ratios  are  substantial  The  mass 
median  diameter  to  collector  diameter  ratio  ranges  from  4.1  X 10-4  to  3.4  X 10'*.  The  effect  of  interception 
would  be  to  substantially  increase  the  values  of  the  mass  median  diameter  efficiency.  However,  with  such  large 
diameter  ratios,  the  collection  efficiency  may  not  equal  the  theoretical  impaction  efficiency  due  to  particle 
bcunce-off  and  re-entrainment  from  the  collector  surface. 

The  Wcibull  distribution  function  is  a simple  mathematical  function  that  adequately  fits  the  data  of 
Landahl  and  Herrmann.  It  may  not  apply  to  all  particle  size  distributions  of  interest  but  as  stated  by  Professor 
Wcibull,  "The  only  practical  way  of  progressing  is  to  choose  a simple  function,  test  it  empirically,  and  stick  to  it 
as  long  as  none  better  has  been  found.”13  The  graphical  technique  for  estimating  the  Wcibull  parameters  as 
presented  by  Nelson15  was  found  to  be  convenient  and  simple  to  use.  and  best  estimates  of  the  three  parameter 
values  can  be  rapidly  found  by  digital  computer.16 

The  computer  program  prepared  for  this  study  and  presented  in  the  appendix  enables  impaction 
efficiency  predictions  to  be  made  over  practically  the  entire  region  of  interest  of  inertia!  impaction  on  cylinders. 
Four  optional  computational  methods  arc  provided  for  the  user  depending  on  (lie  amount  of  input  information 
available.  Option  1 provides  the  impaction  efficiency  of  a cylinder  and  assumes  unit  density  particles  and  an 
ambient  temperature  of  20°C.  The  only  input  requirements  arc  the  mean  windspoed.  cylinder  diameter,  and 
particle  diameter  in  ccntimetcr-gram-sccond  units.  Option  2 makes  use  of  all  the  variables  required  to  compute 
the  K and  <p  parameters  defined  by  equations  3 and  4,  but  the  program  provides  the  necessary  calculations  based  on 
these  mput  values.  Option  3 assumes  that  the  user  has  computed  the  K and  <p  parameter  values  from  equations  3 and 
4 and  desires  only  the  predicted  efficiency.  Thus,  only  the  K and  input  values  arc  required.  Option  4 considers 
the  Wcibull  distnbution  function  and  requires  all  input  data  from  option  2 except  for  the  particle  diameter  for 
which  the  Weibull  parameter  estimates  of  rj.  (3.  and  7 are  substituted.  The  program  computes  the  average 
impaction  efficiency  of  the  distribution,  as  well  as  the  efficiency  based  on  the  mass  median  diameter  of  the 
distribution.  Option  4 has  been  separated  from  the  main  program  so  that  the  reader  could  substitute  his  own 
distribution  function,  if  required,  with  only  a few  minor  program  changes.  A detailed  explanation,  sample  input 
cards,  test  cases  with  sample  outputs,  and  the  complete  program  arc  provided  in  the  appendix  for  the 
convenience  of  the  reader. 

The  approximation  methods  developed  in  this  study  to  calculate  the  impaction  efficiency  of  a 
cylinder  for  any  K and  $ value  have  reduced  the  computational  time  for  the  pomt-by-poinl  trajectory 
calculations  used  to  develop  the  inertial  impaction  theory  from  5.6  seconds  to  13  milliseconds  per  data  point  on 
a digital  computer  comparable  to  a Univac  1 108  computer  without  sacrificing  the  accuracy  of  the  prediction 

V.  CONCLUSIONS. 

1.  Standard  inertial  impaction  efficiency  curves  for  cylinders  developed  from  point-by-point 
trajectory  calculations  have  been  fitted  by  unique  exponential  functions. 

2.  Lagrange’s  interpolation  formula  has  been  adopted  and  applied  to  compute  the  impaction 
efficiency  of  a cylinder  for  any  parameter  values  of  0.13  < K < 300  and  0 < <P  < 10.000. 

3.  The  accuracy  of  the  impaction  efficiency  computations  by  use  of  the  digital  computer 
program  is  withi  1 7<  of  the  point-by-point  efficiency  predictions. 
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4.  The  cumulative  distribution  function  of  the  Weibull  distribution  has  been  incorporated  into  j 

the  digital  computer  program  and  has  been  used  ) fit  the  experimental  particle  size  distributions  of  Landahl  and  1 

Herrmann. 


5.  The  average  impaction  efficiency  for  the  particle  size  distributions  and  the  mass  median 
diameter  impaction  efficiencies  have  been  compared  to  the  experimental  data  obtained  by  Landahl  and  Herrmann 
for  heterogeneous  aerosols. 

6.  Use  of  a mass  median  diameter  to  characterize  a particle  size  distribution  does  not  accurately 
predict  the  experimentally  obtained  impaction  efficiency.  For  large  values  of  the  inertial  parameter,  the  mass 
median  diameter  overestimates  the  efficiency  and,  at  small  inertial  parameter  values,  the  mass  median  diameter 
significantly  underestimates  the  impaction  efficiency. 

7.  Consideration  of  the  particle  size  distribution  in  impaction  efficiency  predictions  results  in 
good  agreement  between  theory  and  experiment. 
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APPENDIX 


DIGITAL  COMPUTER  PROGRAM  FOR  INERTIAL  IMPACTION  EFFICIENCY  COMPUTATIONS 


I.  INTRODUCTION. 

The  computer  program  IMPEF  and  its  subroutines  LAGRNG,WEIBUL,and  POLY  have  been  written  for 
the  Edgewood  Arsenal  Univac  1 108,  time-sharing,  multiprocessor  system. 

IMPEF  has  four  executable  options.  Three  of  these  enable  the  user  to  input  various  particle  and  fluid  flow 
conditions  to  obtain  an  inertial  impaction  efficiency  value.  The  fourth  option  computes  tire  average  impaction 
efficiency  for  a particle  size  distribution  generated  from  the  Weibull  distribution  function  and  provides  the 
impaction  efficiency  for  the  mass  median  diameter  (MMD)  of  the  distribution. 

II.  MATHEMATICAL  BASIS  OF  PROGRAM  IMPEF. 

The  first  step  in  the  development  of  the  approximation  methods  has  been  to  employ  the  point-by-point 
particle  trajectory  method*  to  obtain  the  inertial  impaction  efficiency  curves  for  14  selected  velocity  field  scaling 
parameter  values  chosen  as  standards  (0  = 0;  1 : 2.2;  5;  10;  22;  SO;  100;  200;  500: 1 ,000;  2,210;  5,000;  10,000). 

Polynomial  equations  have  been  developed  that  accurately  fit  the  standard  0 curves.  Each  curve  has  bee. 
sectioned  into  four  parts  representing  four  intervals  of  the  inertial  parameter,  K (0.13-0.22;  0.22-0.5;  0.5-1. 0: 
1.0-64.0).  The  polynomial  coefficients  are  stored  in  the  subroutine  POLY.  If  the  K parameter  lies  between  0 and 
0.13,  the  impaction  efficiency  is  assumed  equal  to  zero  based  on  results  obtained  from  the  point-by-point  trajectory 
model  computations. 

In  order  to  approximate  the  impaction  efficiency  for  any  nonstandard  0 value  over  the  range  1 < 0 < 
10,000,  Lagrange’s  interpolation  formula  has  been  applied  after  assuming  a functional  relationship  exists  among 
the  efficiencies  at  a given  K value.  The  advantage  of  Lagrange’s  interpolation  routine  is  that  equidistant  values  of 
the  independent  variable  are  not  necessary  to  obtain  a reliable  result.  A linear  interpolation  has  been  used  for 
those  cases  where  the  0 value  falls  between  0 and  1. 

Lagrange’s  routine  is  utilized  with  alternating  standard  0 curves  (1;  5;  22;  100;  500;  2.210:  10,000)  to 
compute  efficiencies  for  inertial  parameter  values  in  the  extrapolated  range  of  64  < K < 300. 

The  adaptation  of  Lagrange’s  interpolation  formula  used  in  the  LAGRNG  subroutine  of  the  program 
IMPEF  is  written  as: 


M 


= 


M 

n (In  0 - In  0:) 
E W 

* 5 

n (In  0j  - in  0:) 
j=l 


i = 1 


‘Stucmpflc,  A.  K.  EATR4705.  Impaction  Efficiency  of  Cylindrical  Collectors  in  Laminar  and  Turbulent  Fluid  Flow.  Parti. 
Inertial  Impaction  Theory.  March  1973.  UNCLASSIFIED  Report. 
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where 


= 1 .0/  |exp[exp(Ajyj(!nK)^  + AN.,(lnK)I'M  + ...  + A2(lnK)2  + AjlnK  + AQ)]} 
and  K = inertial  parameter  value  of  interest;  0.13  < K < 300. 

If  0.13  <K  <64, 

= 1;  2.2;  5;  10;  22;  50;  100;  200;  500;  1.000;  2,210;  5,000;  10,000 
i,j=  1,2, ....  13 
M = 1,2, ....  13 
N = 4,  5,  6 or  7 
If  64  < K < 300. 

<>;,  = 1;  5;  22;  100;  500;  2,210;  10.000 

i,  j = 1.  2 7 

M = 7 


A^j's  are  the  coefficients  obtained  from  the  fit  of  the  standard  <>  curves  and  are  given  in  the  subroutine  POLY. 

111.  EXPLANATION  OF  INPUT  FOR  PROGRAM  1MPEF  WITH  EXAMPLES. 

Option  One. 

Computes  the  velocity  field  scaling  parameter.  <?,  the  inertial  parameter,  K,  and  the  inertial  impaction 
efficiency  assuming  an  ambient  temperature  of  20°  C and  unit  density  particles.  Requires  input  of  pat  tide 
diameter,  cylinder  diameter,  and  frcc-stream  wind  velocity. 


format  for  Input  Cards 


Card 

Column 

F<  rmat 

Svmbol 

t ■ - 

Denotation 

Unit 

1 

1 

W 

10P 

Option 

2-80 

79A1 

ID 

Data  identification 

2 

1-10 

FI  0.0 

U 

Frce-strcant  velocity 

cm/sec 

11-20 

FI  0.0 

D 

Collector  diameter 

cm 

21-30 

FI  0.0 

PDM 

Particle  diameter 

microns 

Appendix 
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Sample  Input 


i i i i i * i ii  i i ii  i i ii  iiit  i i ii 


Iramra  npamagroa  rc-iri35XiuO  cjuliJLA±^L.zs^zLria^.r^XT::^x^J'r'^rrzrc  xncgt'L  ».ti.  uti  y 

; 


Sample  Output 


TEST  CASE  FOR  OPTION  ONE 

PARTICLE  DENSITY  = 1. 0000*00  GK/CC  FLUID.  DE  N5IT  Yl  = 1.2047-03  GMyCC 

FLUID  VISCOSITY  = 1.8100-04  POISE  COLLECTOR  DIAMETER  = 1 .4000  CM 

FREE-STREAM  VELOCITY  = 424.650  CM/SEC 

PARTICLE  DIAMETER  = 55.Q0Q  MICRONS 

VELOCITY  FIELD  SCALING  PARAMETER*  PHI  = 42.9023  • 

INERTIAL  PARAMETER*  K = 5.6326 

INERTIAL  IMPACTION  EFFICIENCY  = .7353164  . 


Option  Two. 

Computes  the  velocity  field  scaling  parameter,  <p,  the  inertial  parameter,  K,  and  the  inertial  impaction 
efficiency  for  the  applicable  input  conditions.  Requires  input  of  particle  diameter  and  density,  col’ector  diameter, 
fluid  density  and  viscosity,  and  free-stream  wind  velocity  at  the  applicable  temperature. 


Reproduced  from 
best  available  copy. 
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Format  for  Innut  Card 


S 


Denotation 


2TE5T  CASE  FOR  OPTION  TWO 


hr | i '4  I_n_3f  fr^'ii  V:‘*3  ii^IVii  ii*ii  n~n  i)  ;i  nil  ii'n  II  n U ii'iThTi  »)»  ftjaTTajjii  «"o  « V>  ii  «t_4j  u jiiflgB  g g lilTli  IIM 

*'*  • 1 i~  < J ~ bf  'J  i'2  1 1 ’~S  i>~ ~i*  ~ cl  'n  ;v  joiiV?)!,*  Vo  41  ti  41  aa  aa 


1 1 Ijl  1 I I l Ijljl  1 Ijl  1 ijl  I Ijljl  I Ijl  1 I I 1 ] I 1 1 IJ1  I I t 1 Ijljl  I l|l  I ijl  I I I 1 I I I I I I I I I 1 I I I 1 I I I 1 1 I I I I | 1 1 1 t | 
ife2  2|tt 2 212 2 2 1 2 2 2j2  2 2|2  2 2:2|2  2 2;2  2 2 2 2 2 2 2 22|2  2 2’2  2 2!2  2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 


3 |)  3,3  8 3,3  3 3 3|3  3 3 3 3 3 3 3 3,3,3  3 9:3  3 3 3 3 3 3 3 3 3 3 3 3j3  3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3,3  3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 s . j 
I I I I j 5 H 

14  4 4|4  4 4 4 4 4 4|4  4 4 4 4 4j4  4 4 4 4 4 4j4  4 4j4  4 4j4j4  4 4j4  4 4|4  4 4 4 4 4 4 4 4 4 4 4 4|4  4 4 4|4  4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 | / 

5 5 355  5|5  5 5j^5  5 5j5  5 5j5  5 5|5|3  5 5j5  5 5j5  5 5j5|5  5 5|5  5 sjs  5 5 5 5 5 5 5 5 5 5 5 5j5  5 5 5j5  5 5|5  5 5j5j5  5 5,5  5 5,5  5 5 5 5 5 5 5 S 5 5 5 5 5 ° \ j 

4JL6J  6 6L6J.lii:6 S S.6  5 8 6 6 6 jj  6 6 siS  8 6 6 6 6J:6  6 6;6  6 6.6  6 $1^6  6 6 6 6 6 6_6  6]6  6 6 6J_6  6]6  6 6.S  6 6_il$  6 616  6 |i6  i_6  6J6  6 6 6 6 616  « | \ i 

, ONE  | TWO  1 THREE  ' FOUR  I * FIVE  ‘ " SIX  ' " SEVEN  EIGHT  JJ  3 fc| 

7 7 )|r?  7l7 1 ?j)i| 7 7{7  7 7p 7 7j2jO  IjT T7|7  T7iljf7Tj)  TTjl  7 7|7p  7 7j7  7 7|7 77j7  7 7 7j7  7 7j7  7 717  7 7 7j7  7 7|7  7 717  rTTTTTTTJO  7p  : ^ 

ft  0 ft  ft  ft  ft  Q 9 ft  a a ft  ft  ft  r Q ft  ft  F ftLfi  ft  a,Q  ft  c ft  ft  0 9:9  e ft  0 0 cio  ft  ft  aft  0 ft  : 0 oL  e o,o  0 ■ ft  a 0 o a a do  a e . e * aL  a o'o  o a o a a o a a aL  . 


is  9 99  S sis  9 909  « 9*8  9 99  9 B99 , 9'9  9 9j9  9 sl}s  9 9!9  9 g|9  9 9|9!,  9 9'0  9 9i9  , ^ ^ 5 J1  H 9 I } I g 9I9 

,1  1 a I a j <1,  1 1 '8 ,1 1.  « ’amt  1 i :i  :>  .*  n .v ;i  n *;n  u ujy  a m]i,  a<  *;««  *1  eaa « u[aa  <4  «be|ii  uujaauu.a.  sasa^Jat  u u!*a  e nja>  aa  uj»fti  n »|»  n ?»j„  it  n|u 

rttaiiaiavia 


Sample  Output 


TEST  CASE  FOR  OPTION  TWO 

PARTICLE  DENSITY  : 1. 0430*00  GN/CC  FLUID  DENS  IT  Yl  = 1.2047-03  6M/CC 

FLUID  VISCOSITY  = 1. 8100-04  POISE  COLLECTOR  OIASETER  = 1.2300  CM 

FREE-STREAM  VELOCITY  = 312.500  CM/SEC 

PARTICLE  DIAMETER  = 75.Q00  MICRONS 

VELOCITY  FIELD  SCALING  PARAMETER*  PHI  = 25.9459 

INERTIAL  PARAMETER*  K = 9.3793 

INERTIAL  IMPACTION  EFFICIENCY  = .8221419 


Option  Three. 

Computes  the  inertial  impaction  efficiency  for  a given  set  of  K and  $ input  parameter  values. 

Format  for  Input  Cards 


Card 

Column 

Format 

Symbol 

Denotation 

1 

1 

11 

IOP 

Option 

2-80 

79A1 

ID 

Data  identification 

2 

1-10 

F10.0 

PHI 

Velocity  scaling  parameter 

11-20 

FI  0.0 

RK 

Inertial  parameter 

Sample  Input 


8.4367  23.9025 

i i I ii  rw#i,^Tl  i i it  i i ii  i i ii  i i m i i ii  i i it 


Imam  unagiraiixmi 


i oddi  i ni  a wzsacma  oieb-  u 


)Tni[rrmn7imim 


9 9 9,9  9 9 9 9 919  9 9 9|9  9 919  9 9,919  9 919  9 9,9  9 9 


T CASE  FOR  OPTION  3 

OCITY  FIELD  SCALING  PARAMETER*  PHI  = 
RTIAL  PARAMETER*  K = 23.9025 

RTIAL  IMPACTION  EFFICIENCY  = . 922  36  95 
lNTIAL  FLUIO  FLOW  MAY  NOT  APPLY 


8.4367 


GENERAL  PURPOSE  0 fiELD 


b 6 6]6  6 6L6j6  66166  5 


lrmiTuriTOTi 


I 


TEST  CASE  FOR  OPTION  FOUR 

PARTICLE  QENSITY  S 1.0*JO*00  6M/CC  FLUID  OENSITYi  s 1*20*7-03  GM/CC 
FLUIO  VISCOSITY  = 1*8100-0*  POISE  COLLECTOR  DIAMETER  = .9000  CM 

FREE- STREAM  VELOCITY  r 357*600  CM/SEC 
VELOCITY  FIELO  SCAL2NS  PARAMETER.  PHI  = 22.2678 

VEI3ULL  PARAMETER^  SCALE  = 3*.*356  SHAPE  = 1.7*61  LOCATION  r -1.7375 

AVERAGE  INERTIAL  IMPACTION  EFFICIENCY!  = .*666792 


CUMULATIVE 

MAXIMUM 

INTERVAL 

INERTIAL 

°ER  CENT 

PART  OIAM 

AVERAGE 

INERTIAL 

IMPACTION 

HASS 

(MICRONS) 

(MICRONS) 

PARAMETER 

EFFICIENCY 

2 

2.13* 

2.6*7 

.0178 

.0000000 

3 

3.160 

3.607 

.0331 

.oooonoo 

* 

*.055 

*.*59. 

.0506 

.3000000 

5 

*.86* 

5.238 

.0698 

• 0000090 

6 

5.612 

5.96* 

.0905 

. 3000000 

7 

6.315 

6.6*8 

.112* 

,3000000 

8 

6.982 

7.300 

.1356 

.aonrm 

a 

7.619 

7.925 

.1598 

.3015139 

10 

8.232 

8.528 

.18  50 

.3071181 

11 

8.82* 

9.1 12 

.2112  • 

.0161853 

12 

9.399 

9.679 

.2383 

,0288266 

1 3 

9.959 

10.232 

.266* 

.0*29852 

1* 

IQ .506 

13.773  • 

.2953 

.058  •'635 

1 5 

11.0*1 

11.303  • 

.3250 

.07*8713 

16 

11.565 

11.823  • 

. 3556 

.3916102 

17 

12.081 

12.335 

. 3871 

.108*189 

18 

12.589 

12.8*0 

. *19* 

.1253332 

19 

13.090 

13.337 

.*525 

.1*22780 

20 

13.585  ■ 

13.829 

.*865 

.1585318 

21 

1*  .073 

1 * . 3 1 5 

.5213 

.17*9852 

22 

1* .557 

1 *. 797 

.5570 

.1913652 

23 

15.037 

15.2  7* 

.5935 

.2068813 

y 
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CUMULATIVE 
R£R  CENT 
MASS 

24 

25 

tL  6 

27 

2 s 

2? 

7 

31 

32 

3 * 

34 
25 
3 5 
37 
3 r 

3 9 
49 

4 1 
4 J 
4 2 
4 4 

45 

46 
4 7 
4 8 
4 3 


MAXIMUM 
PART  OTAM 
(MICRONS) 
15.512 

15.964 

IS. 453 

IS. 319 

17.363 

17.34b 

i c . 3 r.s> 

1 o • 7 oS 

1 a . 2 2 i 
19.651 
2a.  138 
2J.SS4 

1 . u 5 1 
.1*5 uo 
2l . 9bo 
: 2 .4  24 

2 : . 3 c 4 
.3  . 34 j 

23.6  27 

44.271 

24 .7  27 
25.235 
25.575 
25.148 

26 .8  <:S 
27.134 


INTERVAL 
AVERAGE 
( MICRONS  ) 

15.748 

16.2  13 

16.6  86 
17.1  5.1 

17.6  14 
16.376 
16.525 
lf .534 
19.452 
19.909 

23 . 3 oo 

23.6  23 
2 1 . 2 o3 
21.737 
22.1 35 

22.6  54 
23.114 
2 3.-76 

24.3  39 
24 . 5 u 4 
24.971 
25.441 

25.812 
26  . 7 87 
26.564 
27.345 


INERTIAL 

PARAMETER 

.6205 

. 65  31 

. 73b2 

.74  84 

. 7s  9: 

.9212 

. 37  43 

. 91  75 

. s=  2i 

1 .3354 

1 . 35  52 

1 .1.-3" 

1.15 

1 . 23  29 

1 . ■’5  2 2 

1. 335= 

1 . 35  92 

1.414' 

1 .4701 
1. 527= 

1 . 5s  6. 7 
1.64fac, 
1.7361 
1 . 77  13 
1 • ^ 3 S'* 

1 . 952? 


INERTIAL 

IMPACTION 

EFFICIENCY 

.2219787 

.2369364 

.2518'* 21 

.266 5775 

. 28P5°n7 

.2348335 

. 3381 317 

. 3 23  78  12 

. 3 

. 3 J 5 ' '5 
. 35= 14=5 
. ? 73  3 *»  11 
. la' 

• 333'.cll 
.4352972 
.41o373c. 
.4271942 
. 4 37  77  13 
.4461161 
. 4 562405 
.4e81c46 
.4773582 
, 4 3 7 3 9 n 4 
.4ys  72  9s 
. : 0= 594s 
.514  9=?4 
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Reproduced  from  g. . ■ 

best  available  copy.  s * 


MAXIMUM 
PART  DIAM 
(MICRONS) 


INTERVAL 
AVERAGE 
(MICRONS  r 


INERTIAL 

PARAMETER 


INERTIAL 

IMPACTION 

EFFICIENCY 


50 

27.587 

27.130 

1.9361 

1.9703 

.5193725  MHO 
.5237306 

51 

28.073 

28.318 

2.0401 

.5  324161 

52 

23.564 

28.811 

2.1117 

.5409*:  54 

53 

29.058 

28.308 

2.1852 

.5493546 

5 <4 

29.558 

29.810 

2.2607 

.5576204 

55 

30.062 

30.317 

2.3382 

.565757^ 

5 6 

30.572 

33.829 

2.4173 

.5737727 

57 

31 .087 

31.347 

2.4996 

.5816702 

56 

31  .60S 

31.871 

2.5842 

.58945=5 

53 

32.135 

32.402 

2.6709 

.5971336 

60 

32.669 

32.940 

2.7603 

.6047093 

61 

33  .210 

33.485 

2.8524 

.6121873 

62 

33.759 

34.038 

2.9474 

.6195723 

63 

34  .316 

34.599 

3.0454 

• 6 268  6 86 

64 

34  .882 

35.170 

3.  1467 

.6  340608 

65 

35.457 

35.749 

3.2513 

.6412132 

6 6 

36.042 

36.339 

3.3595 

.6482703 

67 

36  .637 

36.940 

3.4715 

.6552563 

66 

37.243 

37.553  • 

3.5876 

.6621757 

69 

37.862 

38.177 

3.7073 

.66SC7?e 

7 3 

33.493 

38.815 

3.8323 

.6755322 

71 

39 .138 

39.468 

3.962E 

.6825782 

72 

33.798 

40.1 3S 

4.0980 

.6892757 

73 

40  .473 

40.819 

4.2383 

.6959295 

74 

41.166 

41.521 

4. 36  59 

.7025444 

75 

41 .877 

42.242 

4.5395 

.7091256 

'JMIA.  AT  IV  E 
P£ R CENT 
HASS 

MAXIMUM 
PART  OIAH 
(MICRONS) 

INTERVAL 

AVERAGE 

(MICRONS) 

INERTIAL 
PAR  AM*  TER 

inertial 

IMPACTI0N 

EFFICIENCY 

76 

42.608 

C2.984 

4.7C3'1 

.7 156791 

77 

43.360 

43.748 

u.  os*: 

.7222102 

78 

44.136 

44.537 

5.346? 

.7  28  7''c^4 

79 

44.338 

45.353 

5. ■'323 

.77.62  315 

83 

45.769 

4 6.193 

S. 4t33 

.74173c7 

91 

46.630 

47.378 

5.  -.Sc: 

.7462463 

92 

47.525 

47.992 

5. a595 

.7547723 

9 3 

43.453 

45.34b 

c. . "j  3 5 1 

.7  61  32  40 

5,  4 

4 3 .4  3o 

49.548 

6 . 34  65 

.767  M ?e 

u l 

5 J . 4 1)0 

d1  .303 

□.616: 

.774='?? 

3 6 

51.5  35 

d2.1  13 

6. 903' 

.7812=77 

87 

52  .bo J 

53.2  68 

7.  ?2  «. 

. 7 St"  ~4  7 1 

9? 

13.393 

34.541 

7.5: 77 

. 7 94  54  = j 

53 

55.189 

55.5  86 

7.94  57 

.53197^0 

3-3 

58.564 

i7.341 

3.  7o4^ 

.fc-’M??! 

31 

u .036 

dc . 3 27 

5 . 5 3 2 3 

. 5 1 o 5 ”*  1 1 

U 

39.757 

fa  j . f . 7 7 

3 . 38  6 'J 

.3  24  27  7 = 

93 

61.538 

6 2 .£  77 

3. 11 

.b  3?31 °1 

5“ 

6 3 .6  75 

64.673 

10. 7365 

. 8 4(39  7 14 

3 ~ 

66  .071 

67.495 

11.55  93 

.8 5° 7"  in 

3 \ 

69.313 

73.694 

12.7140 

.5  631 

57 

72  .463 

74.663 

14.2399 

. R 7155(:c 

9 c 

77  .268 

- K -■’T  PRIM5 
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CL*UOE.p£LLEGRlNO.IMP£F 


i 

C 

MAIN  PROGRAM  IHPEF 

i-PEF"-'! 

j 

DIMENSION  P0138I .APG(97).RK SI  971 .El  971  .101791 

I .jrrn'i 

3 

IPRTrO 

iMPEFn-” 

«« 

1 

READ  (5.39.£NQ=29>  IOP.IO 

I «'P£FC“9 

a 

IF  (IOP.LT.l.OR. IGP.GT.9)  50  TO  26 

j uprf .**  j, 

u 

GO  TO  12.9*5.1 11.  IOP 

IHPCFO-'' 

1 

2 

READ  15.351  U.O.POM 

Im Pfr"-7 

o 

RHC=1.C 

i “PEF  3“5 

9 

RH0A:1 .2H97E-3 

IMajPfjno 

1 «> 

RKJ=1.81E-9 

1“1>EFQ1  ~ 

11 

3 

RK:( POM*  POM. 1. JE-81 •! (RH0»U1 /19.0»R«U«0)  1 

I MP£F  91  1 

12 

9HI=(3.0.RH0A.RH3A*U*01/(RHU>RH01 

I“p-:fci  ? 

11 

50  TO  6 

iHPCFCl  ’ 

1** 

9 

READ  (5.36)  RHO. RHO A. R“U. U* O.PDM 

. “P£F  01 9 

13 

GO  TO  3 

I«?cf-ii  = 

IS 

5 

READ  (5.38)  PHI ♦ RX 

I-PEr  71  R 

17 

6 

IF  CIPRT1  7.7.6 

i«P£F  31  7 

1: 

7 

WRITE  16.32)  ID 

IHPCF.C1 

1 J 

GO  TO  9 

1 yP£F  01  '■ 

J ^ 

e 

WRITE  (6.33)  ID 

j IO£fO'  « 

21 

c 

IF  HOP. EC. 3)  GO  TO  10 

i«P£F72  1 

2? 

WRITE  (6.91)  RHC.RHO* .RHU.O.J 

I.-Pff  cr? 

23 

WRITE  (6.92)  POM 

; 7 

£ 9 

1C 

IF  (RK.GT.0.13.ANO.RK.LE.300.)  GO  TO  12 

1 wacro'9 

2a 

IF  (RK.LE.O. O.OR.RK.5T.3DO.)  GO  TO  27 

1 xP'FO:. 

25 

EFF=3.0 

: 

27 

GO  TO  22 

l(.pccT>7 

i 5 

11 

READ  (5.37)  RHO.RHOA.RHU.U.O 

I -U?EF 0?  0 

29 

3HI=(9 .u.RH0A.RH0A*U«0)/ (RHU*RHO) 

: “P£F 029 

3-J 

12 

IF  (PHI. IT. 3.0. OR. PHI. C-T.  10000.0)  GO  TO  28 

I M°£F2x^ 

21 

IF  (PHI. LE. 1.0!  50  TO  13 

IMPEND) 1 

32 

CALL  NUHRTR  (PHI) 

IHPEF03? 

23 

13 

IF  ( IOP. RE. 9)  GO  TO  21 

I M°EFQ2  3 

39 

WRITE  (6.32)  ID 

IMPEF039 

J3 

WRITE  (6.91)  RHO. RHOA. RHU.O.J 

I HP  E F 0 3 5 

2o 

WRITE  (6.93)  PHI 

19DEF036 

37 

CHL0:( (RHO*U)/ (9.»RMU»0) >«1.0£-8 

IMPEFO’7 

3? 

CALL  UEI3UL  (CHLD.RKMMU.PO.APD.RKS) 

IMPEF038 

39 

SUME:0  .0 

I hP£F  03  9 

90 

00  16  1:2.97 

I'MPEF09O 

91 

IF  (RKS<  I) .LE.3.1 3)  GO  TO  19 

IMPEF09 1 

92 

CALL  LAGRMG  ( PHI . RK S! I ) . E < I )) 

I9PEF092 

93 

GO  TO  15 

I HP  EF  99  2 

99 

19 

E(I)=0.0 

1 HPEF09  9 

95 

15 

SUHE=SUME»E(I) 

I MPEF09  5 

96 

16 

CONTINUE 

IHPEF096 

97 

AVGE=SUME/96. 

IMPEF097 

93 

WRITE  (6.99)  A VGE 

I MPEF09A 

99 

WRITE  (6.31) 

I HPEF09  9 

sc 

WRITE  (6.95) 

I MPEF05  n 

SI 

LIN£S=12 

IMPEF051 

52 

DO  20  1:2.97 

IMPEF0S2 

S3 

L INES^lINES*  2 

IKPEF053 

59 

IF  ( I. £0 .50)  GO  TO  17 

IHPEF059 

55 

WRITE  (6.96)  I.PDlD.APOlD.RKSm.ElI) 

IHPEF055 

56 

IF  (LINES.LT. 55)  GO  TO  20 

IMPEF056 
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57 

LINlS-9 

lMPfF9t  7 

sa 

WRITE  16.301 

I HPEFOr  a 

59 

WRITE  16.951 

1 Mprp  3 

C w 

SO  T 3 20 

l “PE" OA  "> 

91 

17 

IF  1 RKMHD.LC .0.131  GO  TO  19 

IMPEFO'  I 

c 2 

CALL  LAGRNG  (PHI.RKHM0.FFF1 

I“OEFG6> 

55 

18 

WRITE  16.971  I.P01 I l.RKHHO.EFF.AOOl 11 .RASCIl .F 

(I  1 

I ,“PEFr,6  3 

t 4 

SO  TO  20 

7 uarf  n;  u 

99 

19 

£FF:0.9 

1 Marc  qr 

CO 

30  TO  18 

i “Rff  f 

6 7 

20 

CONTINUE 

iMPEFr  7 

6" 

1:98 

l MPEF  07  0 

b J 

•RITE  16.961  I.P01I1 

I “Pm: « 

7„ 

I?RT=-1 

1 MO  Zf  ~ 

71 

SO  TO  23 

1 HOFF  07 1 

7c 

21 

CALL  LAGRNG  1 PHI . RK • £ FF 1 

I “Pf F 7’ " 

72 

22 

WRITE  (6.391  PHl.3K.EFF 

i MP£F  IX  : 

?4 

IPRTRl 

I m»FF  07  4 

75 

23 

IF  (PHI-10.1  29.25.25 

I “Of F ”77 " 

7b 

29 

WRITE  (6.90) 

j “orf  -7 

77 

25 

WRITE  (6.311 

1 uJfFn7 7 

7c 

SC  TO  1 

j MOfF  07  = 

79 

26 

WRITE  (6.961  IOP 

l uo.-FV  - 

c w 

00  TO  25 

I “PfF 08  “ 

0 1 

27 

WRITE  (6.991  RK 

i MPtnF  1 

7 3 

30  TO  29 

1 “OfF  OF  ? 

cl 

23 

WRITE  (6.50!  °HI 

1 MPf FO* ’ 

r 4 

29 

CALL  EXIT 

i MPE"  T1  b 

03 

30 

FORMAT  ( 1H1I 

XH?EF09r 

oo 

11 

FORMAT  ( /I 

1 MPf  ? 9C  C 

57 

72 

FORMAT  (1H1.79A11 

I “OFFTA  7 

83 

33 

FORMAT  i 1H  . 79A1 1 

IHP£F  9*8 

59 

39 

FORMAT  ( 11.73*11 

I MPf F 059 

90 

3S 

FORMAT  (3FlC.ni 

lMPfF"'-- 

91 

36 

FORMAT  ( 3E10.0/3F10.0I 

I.MPtFC°l 

32 

37 

FORMAT  ( 3F10.0/2F12.01 

1 MPEF  93  7 

93 

38 

FORMAT  ( ?F 10.3 ) 

I m?EF  C9  7 

99 

39 

FORMAT  ( 9 QH  VELOCITY  rULO  SCALING 

parameter. 

PHI  :. Fii.:i/29H 

INTRI MPEFCq9 

95 

1IIAL  PARAMETER.  K : .F 9.9/32H  INERTIAL  IMPACTION  EFFICIENCY  : .F 9. 7 ) X “PEF 09 5 

96 

90 

FORHAT  ( 35H  POTENTIAL  FLJIO  FLOW  may  NOT  APPLY) 

IMPEF096 

57 

91 

FORHAT  ( 19H  PARTICLE  OENSITY  R.IPE10.9.6H  GM/CC*  9X , 15HFLJID  OEN SI T1 MPE F 09 7 

98 

1Y  =. 1PE10.9.6H  GM/CC/18H  FLUIO  VISCOSITY  :.1PE10.9,6H  P0ISt.9X 

t 

IH?ffO°8 

99 

220HC0LLECT0R  0 1AME TER  :,0PF7.4.3H 

CM/23H  FREE 

STREAM  VELOCITY 

- « 

I MPE  fob  9 

1 00 

30PF9 .3. 7h  CM/SEC1 

IMPEF17P 

101 

92 

FORHAT  (20H  PARTICLE  OIAMETER  = .F8 

. 3*  8H  MICRONS) 

I MPEF  mi 

132 

93 

FORMAT  (90H  VELOCITY  FIELO  SCALING 

PARAMETER. 

PHI  : . F 1 1 .4 1 

1 MPEF 102 

103 

9b 

FORMAT  (90H  AVERAGE  INERTIAL  IMPACTION  EFFICIENCY  =.F9.7) 

IMPEF103 

109 

9 5 

FORMAT  ( 3X.10HCUHULATIVE.6X.7KHAXIMUM.6X.  8HINTERVAL. 1 7X* 

IMP£F199 

105 

18HINERTIAL/9X.8HPER  CENT. 6X .9HP ART 

0IAM.SX.7MAVEPA5E.9X. 

I MPEF  ns 

1 C6 

28HINERTIAL.6X.9HIMPACTI0N/oX.9HMASS.8X.3H(MICR0NS).9X.9H(H1CR0NS) 

.IMPEF 106 

107 

33X. 9 HP ARA ME  TER. 9X.1 OH EFFICIENCY! 

IMPFF 107 

138 

96 

FORMAT  ( 6X.  I 3.  cX  . F 8 . 3 /-3  1 X . F 8.  3»3X. 

F9.9.5X.F9.7) 

IMPEF 108 

109 

97 

FORMAT  (6X.I3. 9X.FB.3.16X.F9.9.5X. 

F9.7.9H  HMD/  31  X » F 8. 3 . 3X  »F  9'.9 

• 

if* 

X 

.IHPEF109 

no 

1 F 9 .71 

IMPEF110 

111 

98 

FORMAT  ( 7H  OPT  ION. I 2. 1 7H  IS  NOT  AVAILABLE) 

X MPEF  11  1 

112 

99 

FORMAT  (29H  INERTIAL  PARAMETER.  K 

:,E19.9.1 JH 

OUT  OF  RANGE) 

I MPEF  11? 

1 13 

50 

FORMAT  C90H  VELOCITY  FIELO  SCALING 

PARAMETER. 

PHI  I.FI9.9.13H 

OUT 

I MPEF m 

1 19 

l OF  RANGE) 

I MPEF  119 

115 

ENO 

I MPEF  1 1 5 
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X AUOE  »PELLEGRIN0.W£I8UL 


SUBROUTINE  MEIBUL  CCHLD.RXHHD.PO.AP&.RKS) 

DIMENSION  POI98).APDC97t*RXSI97» 

READ  <S.2»  ETA. 8ETA. GAMMA 
WRITE  16* 3)  ETA.SETA.GAKNA 
POWER: S« 0/9ET A 
OIFF:£TA-GANNA 

PDt  INO:GAHMA*  lOIFF**!  I ALOGI  1«0/1 1 .0-tFLO  AT  t INC!  /109. 11  It  **POMER  )'• 

IF  liNC.E0»2l  GO  TO  1 

I=INC-I 

Af>om=<POtINCI*POU)  1/2. 

RKS(I):*°01  I»*APO<H*CHLD 
CONTINUE 

RKhSD:P0fSa>«P0«50»*CHlD 

RETURN 

fSSS2t  t 28h°uE I8ULL  PARAMETERS!  SCALE  =.F10.*.9H  SHAPE  z.tU*. 
U2H  LOCATION  :.F9.Rl 
END 


WEIBJLTI 
UEI9ULT? 
UEI8ULOT 
WEI3UL1P 
WEI3ULDS 
WEI3ULTS 
WEI3JL~ 7 
UElBULOfi 
WEI9‘JL',9 

WEI  9UL1  *5 
WE18UL1 1 
WEI  8-JLl  7 
WEI9UL1 T 
WEI3UL1R 
WEI3UL1S 
VEI3UL1S 
WE18UL17 
JE13UL1 f 
W El  B’JLI  9 
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CL  AuDE  »P ELLEGR  INO 

i.LAGRNG 

1 

SUBROUTINE  LAG RhG  t PHI • R<  «£FF  1 

LAGRNGOl 

2 

DIMENSION  ENC131  .ED  11  3.1 . 0 IF  1 1 3 1 . X 1 1 31  .ELD ( 7 1 .ELN(  7 1 . DIFL  1 7 1 . Y { 2 ) . LAGRNG^? 

3 

1RNT<'3) 

LAGRNGT3 

*4 

DATA  X/0.0.. 78 8457 356. 1.60943 721. 

2.30258510. 3. 09 13 4245. 3. 91232301. LA GRNG” 4 

5 

14. 60517019. 5. 29833737.6, 214 60 867. 

G. 93775=26. 7. 70 074779. 8. 517193 20. L«5RNGne 

b 

29.21034u38/E2/ 22 35702 0.1 .-1812693 

.4 7. 269447. 688. -8 63 11. 7591 . LAGRNG°6 

7 

341 628.  3997. -2.3 779. 864  3. 18228. 2921 

.-24  986.1752.4  5785.74  57.-  LAGRN5',7 

a 

489020.3004. 330206.731 .-1864839.07 

• 17170? 87. 1/CLO/10098.27S3.-  LAGRNG"? 

3 

5 1523. 10928.6 10. 97 74 5. -4 79.261 7b. 640. 540964 .- 1 SOI .68156.8921 . 767 1 / .LAGRNGT1 

13 

ERNT/Q. 22.0.5.1.0/ 

LAGPNG1  ■> 

11 

KIN  = 1 

LAGO'iGl  1 

12 

1 

IF  IRK.LT.RNTCKIN1I  50  TO  2 

LA3RNG1? 

IS 

XIN=XIN»1 

LA5RNG1  ’ 

14 

IF  1KIN-31  1.1.2 

l ASRNGl 4 

1-j 

2 

RKLX: ALOS ( RX 1 

LAGRNG1 5 

lb 

IF  (»HI. LE. 1.01  GO  TO  5 

LAG Rfl.il  F 

17 

£FF = 0.G 

LA3RNG1 Y 

18 

INC=1 

L A3RVG1 » 

>3 

IF  (RK.GT.64.0) 

L AGPN31 » 

23 

00  4 n=2« 14, INC 

l AjRNG’ “ 

21 

CALL  °OLY  (N.XIN.RXLX .El 

I AGRN3"' ! 

2’ 

IF  ( INC.E0.2l  GO  TO  3 

L A3RN62  2 

23 

EFF=EFF.EN(N-1 l/EOCN-ll.E 

LAGSNG7  l 

> u 

30  TO  4 

L A5RN 3 ' 4 

25 

3 

l-N/2 

l AGRNG’t 

26 

•FF:EFF.ELN«I1/£lO( tl • £ 

LAGP-ig-f 

27 

4 

CONTINUE 

L ASSN G2 7 

23 

RETURN 

L A5RN5? 

29 

5 

DO  6 N-l  . 2 

L AG=NG,'> 

S3 

CALL  POLY  (N.KIN.RKLX.El 

L AC  PNG  7 “ 

31 

Y(NI=E 

LAGRNG7 1 

32 

6 

CONTINUE 

L A3 PNG 1 7 

S3 

EFFSYC 11  •( 1.-PHI l.Y (2 l.PHI 

LAGRNoM 

34 

RETURN 

L A3PNG3 4 

35 

ENTRY  NU4RTR(PhI1 

LA5RNG15 

36 

phil=alog ( PHI  1 

L AGRNG3  c 

37 

00  7 1=1.13 

LAGRNG37 

38 

0IF(I!=(PHIL-X(I11 

LA3RNG39 

33 

IF  < I.GT.71  GO  TO  7 

LAGRNG33 

43 

0IFLII1=(°HIl-X(2.I-111 

LAGRNG4P 

41 

7 

CONTINUE 

LAGRNG4  1 

42 

00  9 1=1.13 

LAGRNG4? 

43 

EN'JH  = 1 .0 

LAGRNG4  3 

44 

ELNUM: 1.0 

LAGSNG4U 

4 5 

00  8 jsl.13 

LAGRNG4S 

46 

IF  ( I.EO.J)  GO  TO  8 

LAGRNG46 

47 

EN‘JM  = OIF(J).ENUH 

LAGRNG47 

48 

IF  ( I.GT.7.0R.J.GT.7I  GO  TO  8 

LAGRNG48 

49 

ELNUM  = OIFLt  JLELNUM 

LAGRNG49 

50 

8 

CONTINUE 

LAGRNG50 

51 

EN ( I > =ENUN 

LAGRNG5 1 

52 

IF  tI.6T.7l  GO  TO  9 

LAGRNG52 

53 

ELNt I 1 =ELNUM 

LAGRNG5 3 

54 

9 

CONTINUE 

LAGRNG54 

55 

RETURN 

LAGRNG55 

56 

ENO 

LAGRNG56 

• 
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:l  AUDE  .PELLEGR  ino.°oiy 

1 SUBROUTINE  POLY  IN.KIN.X.E) 

2 DIMENSION  NCOF  156) »P< 1 4 • 1 «8J 

* DATA  NCOA/6*  7*  7*8»6*6*6t7<o*7i7i7i6<7t7«8*6»6*7tT»6*7*7«7»6t6»7*?» 

4 lS*C*&i8i&«7i7i8*6«6t6t7i3«4«7<A*Si7*7t8iS*7i7*6tS*8i7*8/ 

5 C »HI  = 0 

6 DATA  (P( I. l.Il . 1 = 1. 6) /-1 2.200402. -35. 459569. -297.1  8542.-462.32274, 

7 1-360 .0659 8.- 111. 552 55/ IP (1.2. J) . J=1 . 7 > /62. 346075 . 355 .9 941 7 . 

i 2 1030. 945 1.14 08.1 269. 1064. 5792.42 1.75102. 64. 7841 330/1 P( 1.3« K).x:l . 

■i  371  /-62. 940926. -1 20.18487. -85. 1649 26. -27. 5 96 7 19. -3. 95 10734,- 

1-  4. 97751282.-3. 91818?7£-2/IP(l»4.L)«L=1.8)/-l.CS36 35 4C-5.1. 93320 78E- 

11  54.-1 .615291SE-3.8.088S181E-3.-2.2809549E-2.4.0923174E-J,- 

12  6. 79783478. -3. 91813S3E-2/ 

13  C PHI  r 1 

14  DATA  (Pt 2. 1. I) .1=1.61/4.984 5675.52. 832639 .21 n.594”‘r. 403.96104. 

1=  1376. 12181. 137. 7212/I3(2. 2. Jl.J  = 1.61/-3.6211216£-2.-  . 17  775906.- 

10  2. S4651 109.-. 75 11 54 14. -1.286 63 31. -.137621 36/1 » (2. 3.  K I .K  = l . 6 )/ - 

17  3 1.7649816. -.283381 6461 .-1.4690678 .-.2281 1208, -.7 7448 736.- 

19  41.1tl9515E-2/(P(2.4.LI.L=1.71/9.3879216£-S.-1.4oSl483£-3. 

j i 59. 629289 2E- 3. - 2. 92444 S S£ - 2. 1.3 17s 18 8E- 'A. -.7 8 75 66 38,-1.  11  1954 5E- 2/ 

::  C = HI  ; 2.2 

< 1 DATA  ( PI  3*  1.  1) . 1 = 1.6) /-26.20a 198.- 2 12. 34232.-668.74366.-1 116.7436. 

22  1 -935 .75446. -29 3.  1 5808/ 1 o I 3.2. J) .J:l. 7) /- 2. 9697$8S. - 1 9. 54 b25S. - 

11  252. 775889. -74. 95 2536. -68. 95 0642, -25.1  35597. -4.  1293537/ IP  I 3. 3. K).k: 

24  31,71  /I .6140825.1  .08 53 372. -.96 742b  S3. -.9300 ’1  28.-. 15264 606.- 

23  4.76202787. 2. 4042651 £-4/1 P( 3.4. L) «L=1. 7)/ 2. 03 184 35C-S.- 5. 9103 702 £- 

io  S4. 5. 59 30 55 2E- 3. - 2. 1 0.7066 3£ -2. 1 .237 59 39 E- 2. -.778542 55. 2. 4 0 63^576-4/ 

27  C »HI  = 5 

23  OATa  (PI4.1. I). 1 = 1, 61/ -2 6. 212424, -214.28974,-701  . 6 1204. -11 48.349.- 

29  1 94  1.3 1809. -307. 88626/ I PI  4,2. J) ,J=1,7)/. 45517389, 2. 6995224. 

30  26-5915579.8.  366522 3. 5. 9148874 .1 . 4 3647 3 1.  . 357 3686/1  ° I 4 , 3.  K ) ,K =1 .7)/ 

31  3 -5. 8 102354.- 14.8 262 74  , -I  3. £2912  3,-5.4  8 03251.-. 8580 77 2.-.  78t>310e9. 

32  41.7  312087E-2/IPI4.4.L) ,L= 1 . 8 ) / - 3 . 61 1 9607E- S. 6.  3913  339E-4.- 

33  54.7997086E-3.2.0067387E-2.-4.7588177E-2.3. 7991 4 1 3E-2 .- .7 77C8 322 . 

34  6.01731202/ 

35  C PHI  = 10 

36  OATA  I PI  5.1.1)  .1  = 1.6) /-SO. 664987. -423. 648 33. -1  4 1 5.  64  74.-2361  .2307  . 

37  1 -1 965. 81 64. -652. 49849/ IP  I 5»2. J).J  = 1. 6 )/-. 062 9683 92. -.19555374. - 

38  2. 240  3795 3.. Oil  24 4 131- -.6091 1288.. 097500309/1  PI  5. 3. K ) . K =1 .7)/ 

39  3 1. 45  97 12  IE-2. 9. 49348899E-2.. 28102331.. 32386  49.. 205  96  467, - 

40  4 .7  6590523.  3. 7057228E-2/I P 1 5 .4.L ) ,L= 1 . 7 ) /9. 66  83056E-S .- 1. 6 339242E- 

4 1 5 3. 1. 1191 8S6E-2.- 3.593  S414E- 2.3. 3364  79E-2.-. 76624,05. 3. 73571 91E- 2/ 

42  C PHI  : 22 

4 3 DATA  I PI  6. 1. 1) .1  = 1.6 )/-76. 950465 ,-651. 69666. -2 206. 645.-3  732. 326  3.- 

44  13153 .5891,-1063. 8574/ I P( 6. 2 . J) . J=l . 7) / 5. 7644622. 37.522245. 

4 5 2100.76219.  14  2.  690  ll  .1 12. 4 3496 . 4 5.959403.  8. 0517  84  2/  IP  (6. 3 .K  ),  K=  1 ,7  ) 

46  3/- 4. 7767558, -12. 444464.-H .500556 , -4. 5 95 0604 ,-.6898708 1.- 

47  4.7M831 659. 6, 6480 331 £-2/ (P( 6. 4, L) ,L=1 . 7 >/ 1 . 71 7T586E-4 . -2 , 6882 966E- 

4 8 5 3.  1.67  804  65E-2.-4. 98471  1 6E -2 . 4 . 9792 63 8 E- 2 . -. 7570784 1 .6 .648065E -2/ 

49  C PHI  = 50 

50  DATA  I PI  7.1.1) . J = 1.6) /-5 9. 93828 2. -501. 29 199. -1675.  1358. -27 93. 8696. 

51  1 -2 325. 6901. -77 1.90325/ IP  17.2. J».J= 1»6)/-. 401 5392 8. -1.7185120. - 

52  22. 6991 33 3. -1. 6 26 66 52. -.89185129,. 21709 08 7/IPI 7. 3. K),K= 1.7)/ 

53  35. 7990554E-2..  31 66300 7.. 67074  199 ,.61650862.. 30 0368 99. -.6 622 564 8. 

54  4 .1 0659907/ 1 PI  7 ,4 ,L) .L  = l. 71 / 1 .8537 19E-4 .-2 .79768S3E- 3. 1 .7095729E- 2. 

55  5 -5 . 069 296E - 2. 5.20 31 50 5E- 2. -.7 3 98 5 798,. 10 65930 3/ 

56  C PHI  = 100 
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57 

58 

59 
53 
51 
62 
53 
64 

55 

to 

57 

59 

59 

7G 

71 

72 

73 

74 

75 
7 o 
77 

76 
79 

£ o 
61 
62 
63 

84 
6d 
86 
67 
6 8 
e 3 
SC 
51 
92 
S3 

94 

95 

96 

97 

98 

99 
100 
101 
102 
103 
134 
1C5 
106 

107 

108 

109 

110 
111 
112 
113 


OAT*  (PC  8. 1. It .1  = 1 .61 /-I 7. 064 731 .-135.07  757.-4  27.92071 .-677. 221 63.P0L TOO' 7 
1 -535 .96207.- 168.  813  51  / (PI 8.2. J!  = 61  /.  34226461 .2.4416096.  POL 7.0058 

26. 3623 809.  7.  94 664  87. 4. 01 34 799.1. 2 230'1 1 1/ <P(8 . 3.K > . K=1 . 61 /-  P0LY0059 

36. 3096724. -9. 421 4242. -•».  666  382  3. -.7  304 11  92.-.  71 5 5094  8..  14  3574  52/  P0LY0960 

4(P(8. 4, L).L=l«81/9. 00 6O305E-5. -1.27339446-3. 6. 5331 65 6E-3»-  POL Y036  i 

51.2515H4E-2.-3-38124S4E-3.  1.7035  729E-2. -.71  366247..  14 15  74  44/  P0LY33*2 

PHI  ; 200  PCLYOO.'* 

OAT A (0(9.  1.11 .1  = 1.61/ -42. 396161.-355.  31  647. -1192. 88 35.-200 3 .967 8 .POL Y^OF 4 
1-1684.9216.-56  6.20499/ (PI  9.2. JI.J  = 1. 71 /I  .06608  94  .8,7591944.  Pf-LYin^  5 

228.50104  4,47.240843.4  2.3  93325. 13.867897.  3. 8 321 726/ IP (9 . 3.K 1.  K=  1 ,7iP0LY0356 
3/-39.46U76 1.-83. 45678 1.-60. 74 8291. -20. 64 1773. -2. 95 97241.-  P0LY33- 7 

4.8094 1 524..  1 9835457/ ( P ( 9, 4.L ) .L  = 1.8)/- 3.4 8B1P57E-4  .5 .4  965 1 92E- 3 .-  POLYOO'  i 

53. 465441 r-2. ,i 1210,716.-. 196 95SB 7.. 16083493.-. 7357441. .19805452/  POL  YOGA 9 

PHI  : 500  POLY3070 

DATA  (P(  10.1,11.1:1  .6) /- 12. 020613. -92.566931  .-284.7568  7.-  P0LY->371 

1 436. 60904.-333.62156. -103. 801 33/ CP( 10.2. JI.J51 .61 /-! .4257360.-  POLYOO” 

28. 0514 717,-18.343394. -19. 614001. -11.35725 5. -1.96313«0/( PI  10. 3.K 1 . P0LY0073 
3K:l. 6 1 /4 .2113941 .7,1224666.4,1793245. 1 .01 218 e4 ,-.56606253.  oCLY397u 

4.27458546/(01  13. 4.L 1 .L=l .7 1/1. 3737&83E-4 726557 86-3.1  .46532051  -POL YO 97 5 

52. -3 .7207S01E-2. 2. 6 3698 55E-2.- .65753359. .274 544/  P0LY30?6 

PHI  : 1000  P0LY3077 

DATA  (P(  11.1. 1). I si. 5 1/1 4. 829 38 8. 97. 320297.240. 339 38 .263 . 37799 • »0LYJ07A 

1109.15 04 4/IPI11. 2, J1.J: 1.81/17. 029397. 137. 42034. 473. 41 88 2.  POL  Y307  9 

2895. 1032 2. 987. 87372. 653. 645. 236. 556 92. 36. 7393«4/ (0(11. 3. Kl.Krl. 71  / POLYOO0'' 
34  7.98  3167.  101. 37 274. 80.S6 71  .2  9.54  84  04,4.661602  4.-.  389930  37.  3CL  Y00° 1 

4 .34078  307/(0(1  1.  4.  L >.L  = 1 . 8 1 / -5  .07064256- 5. 4.  89  48  *SC-«  • -1  . 31 8 IS  1 4f -POL  739=-  7 

53.  1.71 5180 7E -3. -2. 9 5582 58E- 3. -1.2  33 57  8 E- 2. -.6 3 34 6349*. 343 7 63/  °CLYOO°  l 

PHI  r 2210  POL Y30” 4 

OATA  (P(  12.1  .1) . irl  .6)/- 63. 354  1 18.  -536.67989.-  18  -8.279c-  P0LY30te 

1 3078.5398.-2604.681  1 • - 879.  7 1 803/ I o ( 12. 2. J I .Jrl  .71/6.44 395  11  . POLYOO0 6 

241  .998908.  11 2.46  264.1  56.04832. 122.5  5226.49.37317  9.  5.  75  393  36/ (“ C 12 .POLYOO0 7 
33.KI  ,Kil,  71  /- 36. 484  32  7 .-77.S321 74. -62. 802403.-24 .132  723.-  POLYOO00 

44. 32 £1 476. -. 9238608. . 4278 15/10(12.4, LI. L:i.8)/-3. 7 36797 36 -6.  PCLY30' J 

S3.62  3304  5r-5.-2. 3078 772E-4 . 2 . 1 29768 7E- 3. - 1 .0  36  7251 6- 2. 2. 60 19 194 E-  P0LY33  -n 
63. -.62493408. .4278 159 3/  POLYOO0 1 

PHI  : 5000  °0L Y3 09 2 

DATA  (°(  13.1  .1) . 1 = 1 .61 /- 10 3. 5 3582  . -8  ft  3. 6o2nS.- 30 15- IF  79.-  ° 0LY"T>3 

151  39.6  04  8.-4  376.  3562.  - 1 4 87 .8046  / 1 ° ( 1 3. 2.  J 1 •.  jrl  , 7 I / 4. 67e83 15  . PCLY3  03U 

231  . 95072 7.  89. 93 34  64. 1 33.  1 1 754, 109. .29784. 46. 487353. A.  84  35497/ (PI 13.P0LY007S 
33. Kl ,K:i. 71/ 34.1 321 66. 69.52464, S3. 28 1683. 1 9. 34 9472 . 3 . 2204 64 2.-  P0LY3036 

4. 3 86 16475..  5 3065892/1  PI  13. 4. LI .L:  1 .81/ -4 . 4 63 6683 E- 5. 7.  1 149  32  76-4  .-o0LY30Q7 

54.70313886-3.1.7822  1 42E -2. -4. 2365 884E- 2. 4. 44 34 3S2E -2  ,-.63295244  . P0LY3C90 
6.53065886/  PCLY3099 

PHI  : 10000  P0LY3100 

OATA  (P(14,l. II. HI. 51 /l 0.993253. 71. 125908.1 73. 31109. 187. 3633.  P0LY0131 

177.14 7214/IPI14. 2, J1 . J: 1 . 8 1 / 10 .6608 98 . 84 . 79884 4 . 26 6. 49 1 79.  POLYPI "2 

2532.49198.587.36483-384.06109.  137.0  3167,  21.4 3P1 96/ (PI  14. 3. K) .K  = l . P0LY0103 
371 /- 16. 2 3 7284, -3 2. 2 69 9. -24. 0022 93. -8. 3 48 9724.-1. 3870594.-  P0LY31n4 

4.69722816. . 61335532/IPI 14.4.L1 .L=l .8! /-I .0977396F-4, 1.60599676- 3.-°OLY3ins 

59. 72 71 265 6- 3« 3. 0 50 9 84E - 2. -5.4 851 7 38E- 2 .4  ,470  41  HE- 2. -.  60 594 7 81.  P0LY0106 

6.61335528/  P0LY3137 

NS=4»( N- 1 l.KIN  ° OL  YO 1 ^ 8 

HrNCOF (NS  I POLYOln9 

RLLRE=P< N.KIN.  1)  P0LY0110 

00  1 1:2. H ° OL  YO 1 1 1 

RLLRE:RLLRE»X*P(N.KIN.I)  P0LY0112 

CONTINUE  POL  YO 1 1 3 


Appendix 


37 


°0LYD114 
B01Y011 
POL  YD  11 


135KBT  PRINT! 


11*  E:l . /( EXB ( EXP( RLIRE  1 ) I 

115  RETURN 

116  END 


J 'JQ 

Apptndix  I J 

01 234567890123*567890123*5678901 23*567890123*567890123*5678901 2345678901234567890123*5678901 2 


7\  *n 


